Quantum limit of laser cooling in dispersively- and dissipatively-coupled optomechanical systems 
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Mechanical oscillators can be cooled by coupling them to an optical or microwave cavity. Going beyond the 
standard quantum noise approach we find an analytic expression for the steady-state phonon number in systems 
where the position of the mechanical oscillator modulates the cavity frequency as well as the cavity line width. 
We trace the origin for the quantum limit of cooling to fluctuations in the optical force both at and away from 
the mechanical frequency. Finally, we calculate the minimal phonon number for the different types of coupling. 
Our study elucidates how to beneficially combine dispersive and dissipative optomechanical coupling. 
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Introduction. Optomechanics is an area of research that is 
concerned with systems in which the position of a mechanical 
oscillator modulates the properties of an optical or microwave 
mode [1-5]. Apart from fundamental questions, e.g. about the 
decoherence of increasingly macroscopic objects [6], these 
systems have promising applications in information and quan- 
tum science, e.g. transducers in quantum hybrid systems [7]. 

To enable these applications reducing the thermal motion of 
the mechanical oscillator has been a focus of intense research. 
Adapting laser cooling techniques from atomic physics [8, 9] 
experiments have observed the quantum ground-state [10-12] 
as well as an asymmetry in the mechanical sidebands [13, 14]. 

In Ref. [15] an optomechanical system has been introduced 
where the cavity line width depends parametrically on the po- 
sition of a mechanical oscillator A strikingly new feature is a 
Fano line shape in the force spectrum that is a consequence of 
quantum noise interference. In a recent study we have shown 
that this form of force spectrum features two cooling and two 
instability regions [16]. Following the proposal in [17] an in- 
terferometer setup has recently investigated these effects [18]. 

For an optimal detuning between laser and cavity frequency 
the zero of the Fano line shape coincides with the mechanical 
frequency. As the quantum noise approach [19] estimates the 
fluctuations of the optical force with the noise spectrum at the 
mechanical frequency, it predicts the unphysical result that the 
phonon number goes to zero as the coupling is increased, i.e. it 
does not give a quantum limit of cooling. 

In this paper we go beyond this level of approximation and 
derive an analytic expression for the steady-state phonon num- 
ber which takes into account the noise in the force spectrum 
at all frequencies. We show that noise away from the mechan- 
ical frequency can become the limiting process for cooling. 
In this case, the steady-state phonon number depends on the 
coupling in a qualitatively different way featuring a minimum 
at finite coupling. We give explicit expressions for the min- 
imal phonon number for so-called purely dissipative as well 
as dissipative and dispersive coupling. While the Fano line 
shape due to dissipative coupling leads to a vanishing ampli- 
fication rate, additional dispersive coupling can increase the 
cooling rate and thus further lower the phonon number Our 
study provides the physical limit of cooling and shows how to 
exploit the presence of these two kinds of optomechanics. 



Model. We investigate an optomechanical system, where 
the resonance frequency Uc of a cavity and its line width k are 
both modulated by the displacement of a mechanical oscillator 
with resonance frequency lom- These two types of coupling 
between optical and mechanical degrees of freedom will be 
referred to as dispersive and dissipative coupling, respectively. 

The Hamiltonian {fi = 1) of such a system is given by 
J^ — (jJcO^a + (jJmc^ c + H^ + H^ + i/jnt, where a^{a) are 
bosonic creation (annihilation) operators of the cavity mode, 
&{c) are bosonic creation (annihilation) operators of the me- 
chanical mode, and H,^ and H.y describe driving and damping 
of the cavity and the mechanical oscillator, respectively. The 
optomechanical coupling is given by [15], 
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where An = — "j^^ ^o is the dispersive and Bn — '^J^^' xq 
the dissipative coupling strength [20]. Here, x — xq{c + c^) 
denotes the displacement of the mechanical oscillator, xq = 
{2mujM)^^^^ is the size of the zero-point fluctuations and m 
the mass of the mechanical oscillator, ¥ ibq) are bosonic cre- 
ation (annihilation) operators of the optical bath coupled to 
the cavity, and p is the density of states of the optical bath that 
we treat as a constant for the relevant frequencies. 

We work in a frame rotating with the drive frequency LUci, 
using a — {a + (i)e"'"''*, where a denotes the mean cav- 
ity amplitude and d the fluctuations around this mean value. 
Then, using the input-output formalism [19] adapted for dissi- 
patively coupled systems, we write down the linearized equa- 
tions of motion [15, 17, 21] 
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with Fxo ~ AHid*d + i^d*^din + i^a*{iA + f )d + h.c. 
Here, A = cj^ — Wc is the detuning between the drive fre- 
quency ojd and the cavity resonance ojc, and din (Qn) describes 
the fluctuations of the optical (mechanical) input mode. We 
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FIG. 1. (Color online) (a) The mechanical spectrum 5cc (i^)i^a/ ob- 
tained from the exact solution (solid black) [16] and the approxi- 
mate expression (5) (orange dotted) for A = 0, B\a\ — 0.2, and 
A — cjm/S. The green dot-dashed (blue dashed) line shows the first 
(second) term in Eq. (5). The thin red solid line is the weak-coupling 
force spectrum Sff (i^)a^o/i^A/ (4). (b) The optical output spectrum 
Sali^) for i = 0, B\a\ = 0.01, A = ujm/2 (solid black), i = 0, 
B\a\ = 0.2, A = ojm/2 (dashed blue), and A\d\ = 0.2, B = 0, 
A — —ujm (green dot-dashed). Other parameters are lom hi = 10^, 
"-th ~ 100, andajAf/«: = 5. 



assume Markovian baths, where the mechanical bath has a fi- 
nite temperature T and thus a thermal phonon number rith = 
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where ks denotes Boltzmann's constant. 



i.e. (cf„(a;)cin(w')) = 27r(5(a; + w')nth and (cin(w)4i('^')> = 
27r(5(a; + a;')(nth + 1), while the optical bath is assumed to be 
at zero temperature, i.e. {^d\^{uS)d\Jy'jj'y) = 27rS{uj + lu'). 

Noise contributions to the mechanical spectrum. For weak 
coupling we can use the quantum noise approach to derive 
transition rates for the mechanical oscillator with Fermi's 
Golden Rule [19]. In our case, the force F leads to transi- 
tions between states with n and nil phonons. The rates 



are given by T„ 
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with the amplification rate F^ — XqSpf{—ujm) and cooling 
rate F^ = XqSff{^m)- They are obtained from the weak- 
coupling force spectrum Sff{<^) = J dt e'''^* {F^ (t) F) evalu- 
ated in the absence of coupling. Here, it is given by [15] 
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As discussed in Ref. [15], for the optimal detuning A = 
iopt — 'jJmI'2 + kA/B a Fano interference leads to a vanish- 



ing amplification rate F^ = XQSFpi—i^M) = 0. The steady- 
state mean phonon number n within the quantum noise ap- 
proach is given by n = (7nth + r-t-)/(7 + F^ ^ T^)- If how- 
ever F^ = and F^ ^ 0, the mean phonon number n goes 
to zero in the limit of large coupling strength, i.e. there is no 
quantum limit of cooling at this level of approximation. 

In the following, we go beyond this standard quantum noise 
approach and take the complete force spectrum SFpii^) into 
account, i.e. the noise at all frequencies. To do so, we solve 
Eq. (3) in the Fourier domain for d{uj) and insert the result 
into the equation of motion of the mechanical oscillator (2). 
Neglecting correlations between the optical field and mechan- 
ical bath, e.g. ((i(aj)cin(w')) — 0, we obtain an approximation 
for the mechanical spectrum Scc{i^) = ^ dt e'^'^*' {c^ {t)c) as 

Scc[i^) = |xa/(-w)P [771th + xISff{^)\ , (5) 

where xm(w) = [7/2— ^(aj—WA/)]^^ is the effective mechan- 
ical response function, and where we have taken into account 
the optically-induced damping 7 = 7 + F^ — F^, but neglected 
the optically-induced frequency shift, i.e. the optical spring. 

In Fig. 1 (a) we plot the mechanical spectrum Scc{i^) for 
purely dissipative coupling A = and detuning A = Aopt. 
Within the quantum noise approach the mechanical spectrum 
has a peak at the mechanical frequency uj = —ujm describ- 
ing the response of the mechanical oscillator to thermal fluc- 
tuations "frith and optical force fluctuations at the mechanical 



frequency, i.e. Scc{uj) = \xm{-l^)\ [inth- 
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In the resolved-sideband regime ojm ^ n the approximation 
(5) features an additional peak due to optical force fluctua- 
tions \xm{^^)\'^XqSff{^) missed by the quantum noise ap- 
proach. For K ^ 7 it is centered at a; = —lum/2 with a zero 
at LU = —COM due to quantum noise interference [15]. In Fig. 1 

(a) we plot the mechanical spectrum iS'cc(w) obtained from the 
exact solution to Eqs. (2) and (3) as given in Ref [16]. The 
agreement with the approximate expression (5) is excellent. 

The reason for the failure of the quantum noise approach 
can be understood by looking more closely at the force spec- 
trum SFpii^) also in Fig. 1 (a). Within the quantum noise ap- 
proach we approximate the optical force fluctuations by evalu- 
ating the force spectrum Sff{^) at the mechanical frequency. 
However, this is only justified if Sff{^) varies slowly around 
zLljm on a scale of 7. For A = Q and A = Aopt this is clearly 
not the case, and the quantum noise approach fails. 

The two contributions to the mechanical spectrum Scc{^) 
in Eq. (5) can also be detected in the optical output spectrum 
5™'(a;) = J dt e^'^* {dl^i{t)dout) where the input-output rela- 
tion is din — ^out — —\/Kd — y/KBdx/2xo [16]. In Fig. 1 

(b) we show the optical output spectrum S^'^{uj) for detun- 
ing A = ujmI'2. It features a dominant peak of width 7 at 
bj = — wa/, a sharp dip to zero at w = ii^A/^ and two smaller 
peaks of width k at w = ±wa//2. The one at w = +a;A//2 ex- 
ists at small coupling, while the one at a; = —ujm/2 appears 
at larger coupling strength. In the limit of strong coupling the 
two peaks have equal weight. This differs significantly from 
the optical output spectrum S°i^{ll:) for purely dispersive cou- 
pling B — Q and detuning A — ~ujm which features the 



well-known mechanical sidebands, i.e. two peaks of width 7 
at cli = iwM, symmetric around the drive frequency. The 
different optical output spectra S'2^'(aj) are a signature of the 
different cooling processes for these two kinds of couplings. 

Improved expression for the mean phonon number Given 
the improved approximation for the mechanical spectrum 
5cc(i^), Eq. (5), we obtain an expression for the mean phonon 
number n by integrating over all frequencies w, i.e. n — 
(c'c) = J /S'cc(w)(ia;/(27r) where we have the analytic result 
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The first term in Eq. (6) accounts for the thermal fluctuations 
due to the mechanical bath with thermal phonon number rith 
reduced by the optically-induced damping. In the limit k ^ 7 
the second term simplifies to XqSff{^'^m)/j, i-C- together 
with the first term it gives the standard quantum noise result. 
Note, that the second term vanishes at the detuning A — Aopt 
and thus it does not provide a quantum limit of cooling. The 
third term in Eq. (6) goes beyond the quantum noise approach 
and is non-zero at the optimal detuning A = Aopt, thus it leads 
to a quantum limit of cooling. For purely dissipative coupling 
A = 0, Eq. (6) coincides with an expression in Ref. [17]. 

For purely dissipative coupling A = at the detuning A — 
uj]\f/2 and in the limit k 3> 7, Eq. (6) simplifies to 
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We see that the phonon number n has a qualitatively different 
dependence on the coupling strength as compared to purely 
dispersive coupling B = Q discussed above. It is surprising 
that the second term in Eq. (7) depends only on the coupling 
strength B\a\ and not on the sideband parameter ijJm / k. 

In Fig. 2 (a) we plot the mean phonon number n as a func- 
tion of the coupling strength B\a\ for purely dissipative cou- 
pling A = {). Within the quantum noise approach the phonon 
number n approaches zero in the limit B\a\ — > 00. For small 
coupling the approximation (6) agrees well with the quantum 
noise approach but, in contrast to the quantum noise approach, 
it features a minimum at finite coupling which defines a finite 
minimal phonon number rimin. Following e.g. Ref. [22] an ex- 
act, but cumbersome and not very illuminating expression for 
the phonon number n can be derived. It agrees very well with 
our approximate result (6) and deviates from it only for large 
coupling strengths close to the unstable region. 

In Fig. 2 (b) we plot the mean phonon number n given by 
Eq. (6) as a function of detuning A. We find that there are two 
cooling regions, close to A = —ojm and around A = lomI% 
where the phonon number is significantly reduced relative to 
the thermal phonon number nth, and two unstable regions 
[16]. The contributions of the second and third term in Eq. (6) 
are shown separately. Notably, the second term in Eq. (6) is 
the limiting term for the phonon number n at most detunings 
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FIG. 2. (Color online) (a) Mean phonon number n obtained from 
the quantum noise approach (green dot-dashed), from Eq. (6) (blue 
dashed), and the exact solution (solid black) for A = 0, u)m /k — 3, 
A — tJAf /2, nth = 100, and coAi/'y = 3 • 10^. (b) Mean phonon 
number n from Eq. (6) as a function of detuning A (solid black). The 
green dot-dashed (red dashed) line shows the second (third) term of 
Eq. (6). Parameters are A = 0, lum/i^ ~ 3, B\a\ = 0.2, nth = 100, 
and um M = 10^. In (c) we show a close up of (b) at A = lom I'i 
and include the quantum noise result (blue dashed) and the exact 
solution (orange dotted). Hatched areas indicate unstable regions. 



A, except close to A = Aopt where it goes to zero. In this case 
the third term in Eq. (6) becomes important and provides the 
quantum limit of cooling, whereas it can be safely neglected 
elsewhere. Fig. 2 (c) is focusing on A w Aopt. It shows that 
our approximate expression (6) leads to a much better agree- 
ment with the exact solution than the quantum noise approach. 
In the following we will focus the case A = Aopt where de- 
viations from the quantum noise result are most significant. 

Minimal mean phonon number For purely dispersive cou- 
pling S = we obtain in the limit k ^ 7 the well-known 



result n = 7nth/7 + xISff{-^m)Ii = (7"-th + r^)/(7 + 
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FIG. 3. (Color online) Minimal phonon number nmin as a function of 
the sideband parameter ojm/k obtained from Eq. (8) (black dashed) 
and the exact result (black solid) for yl = and A — LJi\//2; the 
thin solid black line indicates the limit of rimin for ujm/k — >■ cxd. 
The dot-dashed green line gives the quantum noise result for purely 
dispersive coupling (at A = —ujm, A — > oo), the green solid line 
shows the exact result, optimized for the dispersive coupling strength 
A\a\ (B — 0). The red dotted line shows the minimal phonon num- 
ber for ideally mixed coupling at A = cjai/2 + kA/B, Eq. (9). 
Other parameters are nth = 100 and k/j — IC. 



r^ — r^) from Eq. (6). This leads to a minimal phonon num- 
ber n„in = r^/(r; - r^) = kV(16wm) for A « -ujm and 
A —i' oo. In this case ground-state cooling is only possible in 
the resolved-sideband limit ojm ^ k [8, 9]. 

In the case of purely dissipative coupling A = 0, Eq. (7) 
features a minimal phonon number 
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In contrast to the case A = 0, i.e. Eq. (8), the minimal phonon 
number nmin in Eq- (9) vanishes in the limit ojm ^ k. 

Note that this requires dispersive coupling A and dissipative 
coupling B to have opposite signs. Whereas for either purely 
dispersive or purely dissipative coupling the sign of the cou- 
pling strength does not enter, the relative sign matters in pres- 
ence of both couplings. According to our definitions of A and 
B, opposite signs imply that the cavity resonance frequency 
ujc and the cavity line width k both increase (or decrease) with 
increasing displacement of the mechanical oscillator x. 

In Fig. 3 we plot the minimal phonon number «„;„ for dis- 
persive, dissipative, and mixed coupling as a function of the 
sideband parameter um/k. We see that in the case of purely 
dissipative coupling A = 0, the minimal phonon number «,„;„ 
approaches a finite value for ojj^j ^ k. In contrast, nn,in for 
purely dispersive coupling B ^ goes to zero in the limit 
ujm /k -^ oo. We also show the exact solution and find that 
our approximations capture the qualitative behavior very well. 
Fig. 3 also shows that ideally mixed coupling can overcome 
the cooling limit for systems with purely dissipative coupling. 

In closing we note that opposite signs of dispersive and dis- 
sipative coupling will also allow for cooling on resonance. If 
Ak = —Blum/'^, the optimal detuning simplifies to Aopt = 
ujm /2 + kA/B ~ 0. This would be particularly beneficial 
in the resolved-sideband limit ujm 3> k where off-resonance 
most of the input power is reflected off the cavity. 
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9)^-7(AtV^M + 9)/(4«). 
Remarkably, in the good-cavity limit ujm ^ k, n^i„ becomes 
independent of the sideband parameter lom/i^ and approaches 
a finite value, rimm = ^/^h^^hhJi^^ ~ 97/(16k). 

In systems featuring both dispersive and dissipative cou- 
pling, i.e. A y^ and i? 7^ 0, at A = Aopt, and in the limit 
K ;:^ 7, Eq. (6) again simplifies to Eq. (7). Remarkably, the 
second term of Eq. (6) is independent of the dispersive cou- 
pling strength A, and the minimal phonon number rimin only 
depends on the coupling strength A via the damping rate 7. 

Minimizing Eq. (7) over the coupling strength A leads to 
Ak — — 3i3wA//2. In this case, the optimal dissipative de- 
tuning A — ujai/2 + kA/B coincides with the optimal de- 
tuning for dispersive cooling, A = —lom- Within the quan- 
tum noise approach this means to take advantage of a vanish- 
ing amplification rate F^ = at a detuning where dispersive 
coupling can increase the cooling rate F^ and thus 7. The 
minimal phonon number for this optimal mixed coupling is at 
B^jap = A/nth7(«:Vwif)/«^ - 7('«Vwm)/(4k), and is 



n^th K^ 
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APPENDIX ON DIFFERENT DISSIPATIVE COUPLINGS 

Dissipative coupling arises if the cavity line width k is mod- 
ulated by the mechanical displacement x, i.e. k ~ k{x). So 
far, we have discussed the maximally overcoupled case where 
there is only a single loss channel for photons, so all losses are 
due to the channel where the coherent drive enters. 

In general, several loss channels for photons can be present 
in an optomechanical system. We write the total cavity damp- 
ing K = Koxt + '^o ^s the sum of an external loss rate K^xt 
associated with the channel where the drive enters and kq that 
contains the losses through all other channels. Including the 
additional loss channel leads to an additional input mode a^j^ 
distinct from the input mode af^^ associated with the drive, 
i.e. (a?„) = but (a°^') — a-m- It is this difference that finally 
results in different expressions for the force F below. 

It is crucial to distinguish between the case 



K = Kext{x) +Ko, 



(10) 



where the cavity line width associated with the drive port is 
modulated by the displacement, and the case 



K = Koxt + Koix), 



(11) 



where the internal losses depend on the displacement. 

Case 1. The first case (10) is a generalization of the treat- 
ment in the main text where we used n = Kcxt{x), i.e. kq = 0- 



If Kq 7^ 0, i?Kcxt — '^°dx ^0 leads to the force 
Fxo = Ak (a*d + aS) + iB^/l^t a*d\^^ - (C*)^a 
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uA + ^ j d- ia uA- ^ 



d) 



, (12) 



where d°^^, (df,^')^ are the fluctuations of the input mode, and 
the second line of Eq. (12) is due to the coherent drive entering 
through the same port. This leads to the force spectrum [23] 
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2Ko[(A-2iK/S)2 + «;2/4] 
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(13) 



As discussed in Ref.[23], the optimal detuning Aopt no longer 
leads to an exact zero of the force spectrum Spp (w) due to the 
second term of Eq. (13). Depending on the ratio of Kcxt and 
Kq the quantum noise interference becomes less perfect and 
ultimately, if kq ^ Kcxt, the force spectrum is a Lorentzian. 

Case 2. In the second case (11), where k — Kcxt + ko{x), 
the force differs significantly from Eq. (12). Since the input 
mode associated with ko{x) describes only fluctuations of a 
zero temperature bath, the force is given by 



Fxo = Ak (d*d + aJ^) +^S^/^ [a*< ~ (<)^'a 



(14) 



Notably, a term corresponding to the second line of Eq. (12) is 
missing, since this loss channel is not associated with a drive. 
Thus, for purely dissipative coupling A = 0, there is only one 
contribution to the force and no quantum noise interference. 
The force spectrum simplifies to Sff{^) — KoB'^\a\'^ /{Axq), 
i.e. it becomes completely flat as a function of frequency. 

In presence of both types of coupling, the dispersive cou- 
pling term provides a cavity-mediated force, whereas dissipa- 
tive coupling leads to a force directly proportional to the op- 
tical bath mode rf? , and these two contributions can interfere. 



Using Bkq 



Sppi^) 



dKa{x) 
dx 



xq, the force spectrum is given by 



S2|a|2 Ko[(a; + A - 2Ati/Bf + k^/4] 



B 



4x1 



(k/2)2 + (CJ + A)2 

Kext(2iK/i3)2 



4x2 («;/2)2 + (cj + A)2 



(15) 



i.e. it has a Fano line shape that, in contrast to the first case, be- 
comes a Lorentzian if Kc^t ^ kq. Also in contrast to the first 
case, no perfect destructive interference is possible, i.e. the 
force spectrum Sff{^) has no exact zero. The optimal detun- 
ing from the first case has no longer a special meaning in this 
context and does not lead to a vanishing term of Sff{~'^m)- 
Finally note, that the relative sign of dispersive and dissipa- 
tive coupling, i.e. of A and B, becomes relevant again. 



